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Abstract 

We  present  results  of  experiment  with  two  distinct 
extremum- seeking  adaptive  algorithms  for  control  of  com- 
bustion instability  suitable  for  reduction  of  acoustic  pressure 
oscillations  in  gas  turbine  over  large  range  of  operating  con- 
ditions. The  algorithms  consists  of  a frequency  tracking  Ex- 
tended Kalman  Filter  to  determine  the  in-phase  component, 
the  quadrature  component,  and  the  magnitude  of  the  acous- 
tic mode  of  interest,  and  a phase  shifting  controller  with 
the  controller  phase  tuned  using  an  extremum-seeking  algo- 
rithms. The  algorithms  are  also  applicable  for  control  of 
oscillations  of  systems  whose  oscillation  frequency  and  op- 
timal control  phase  shift  depends  on  operating  conditions, 
and  which  are  driven  by  strong  broad-band  disturbance.  The 
algorithms  have  been  tested  in  combustion  experiments  in- 
volving full-scale  engine  hardware  and  during  simulated  fast 
engine  transients. 

1 Introduction 

Emphasis  on  reducing  the  levels  of  pollutants  created  by 
industrial  gas  turbine  combustors  has  led  to  the  develop- 
ment of  lean,  premixed  combustor  designs,  as  premixing 
large  amounts  of  air  with  the  fuel  prior  to  its  injection  in- 
to the  combustor  greatly  reduces  peak  temperatures  within 
the  combustor  and  leads  to  lower  NOx  emissions.  However, 
premixed  combustors  are  often  susceptible  to  thermoacous- 
tic combustion  instabilities,  which  can  lead  to  large  pressure 
oscillations  in  the  combustor.  These  pressure  oscillations  re- 
sult in  increased  noise  and  decreased  durability.  The  goal  of 


using  active  combustion  instability  control  on  a gas  turbine 
engine  is  to  keep  pressure  oscillations  at  an  acceptable  level 
over  a large  range  of  operating  conditions. 

Experiments  and  model-based  analysis  determined  that 
pressure  measurement  and  a simple  phase-shifting  algorithm 
with  an  appropriately  chosen  control  phase  is  sufficient  for 
suppression  of  oscillations,  given  enough  control  authority 
([11,7,  14,  6,  8,  3,  4]).  Model  based  analysis  determined  that 
minimum  information  needed  to  calculate  the  best  control 
phase  requires  estimation  of  parameters,  including  transport 
delay  (or  at  least  the  corresponding  phase  shift),  that  are 
hard  to  obtain  from  pressure  measurements  alone  [8,  3,  4], 
Even  if  this  difficulty  could  be  circumvented,  the  sensitivity 
to  modeling  errors  was  likely  to  be  high. 

Need  for  developing  an  algorithm  that  would  allow  finding 
the  best  phase  with  minimum  amount  of  a priori  information 
that  would  work  over  large  range  of  operating  conditions  and 
with  minimum  model  assumptions  seemed  apparent.  The 
operating  conditions  include  fast  engine  acceleration  and  de- 
celeration transients. 

An  experiment  with  active  control  of  pressure  oscillations 
in  an  industrial  combustor  has  been  conducted  on  179MW 
and  230MW  gas  turbines  by  Siemens  [16].  The  pilot  fuel  has 
been  modulated  using  a phase-shifting  mechanism.  The  best 
control  algorithm  parameters  were  found  manually  at  each 
operating  condition.  However,  problems  developed  when  the 
conditions  were  changing  suddenly.  Here  is  a quote  from  the 
paper  by  Siemens  engineers  and  they  coworkers  describing 
experiments  with  Active  Instability  Control  on  170MW  en- 
gine: 

” Further  investigations  of  the  effects  of  ATC  showed  de- 
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pendence  of  the  required  control  parameters  (gain  and  phase 
shift)  on  the  operating  conditions  of  the  gas  turbine  (e.g., 
power  level,  pilot  gas  fraction).  AIC  operating  parameters 
wich  result  in  a clear  reduction  of  oscillations  for  certain  gas 
turbine  operating  conditions  can  cause  a significant  loss  in 
performance  of  the  control  system  during  other  gas  turbine 
operating  conditions.  In  the  worst  case,  the  heat  release  rate 
modulation  induced  by  the  controller  can  be  in  phase  with 
the  one  resulting  from  the  self-excitation,  which  results  in 
a positive  feedback  and  further  amplification  of  the  unde- 
sired oscillations.  Therefore,  the  control  parameters  must  be 
adjusted  to  the  operating  conditions  quickly.  Since  the  con- 
trol parameters  were  set  manually,  difficulties  arose  when  the 
operating  conditions  were  subject  to  sudden  changes.  There- 
fore, an  adaptive  control  algorithm  is  being  developed,  which 
seeks  the  optimum  settings  for  gain  and  phase  shift.”  [16]. 

The  lack  of  a control  algorithm  that  would  allow  for  the 
search  of  control  parameters  automatically  and  would  not 
magnify  the  pressure  oscillations  in  the  process  of  tuning  was 
perceived  by  Siemens  engineers  to  be  the  major  obstacle  in 
implementing  Active  Instability  Control  to  suppress  pressure 
oscillation  on  an  industrial  gas  turbine. 

The  algorithm  presented  in  this  paper  allows  fast  auto- 
matic tuning  of  control  parameters  and  (under  reasonable 
assumptions)  can  be  guaranteed  not  to  amplify  the  oscilla- 
tions. Two  extremum-seeking  algorithms  [17]  were  selected 
for  tuning  of  the  control  phase.  There  was  a technical  chal- 
lenge in  finding  a phase  shifting  mechanism  that  works  for 
pressure  oscillations  varying  over  large  frequency  range  and 
reliable  pressure  magnitude  detection  mechanism.  This  chal- 
lenge has  been  overcome  with  a frequency  tracking  observer 
algorithm  (based  on  Extended  Kalman  Filter  described  in 
[12])  that  allows  for  fast  and  reliable  estimation  of  the  in- 
phase  and  quadrature  components  of  the  bulk  pressure  mode 
over  large  range  of  frequencies. 

Performance  specifications  for  extremum-seeking  algo- 
rithms have  been  defined  for  algorithm  initialization  tran- 
sients and  engine  acceleration  transients.  When  initialized 
with  a phase  corresponding  to  amplification  of  oscillations, 
the  algorithms  should  quickly  produce  and  maintain  phases 
corresponding  to  suppression  of  the  oscillations.  In  the  en- 
gine acceleration  transients  the  algorithms  should  be  able  to 
suppress  oscillations  relative  to  uncontrolled  levels. 

The  algorithms  were  tested  on  a single  nozzle  combustion 
rig  at  United  Technologies  Reserach  Center  in  August  1998. 
More  details  on  frequency  tracking  and  classical  extremum- 
seeking  algorithm  are  included  in  [5].  The  paper  [18]  pro- 
vides more  details  on  the  triangular  search  extremum-seeking 
algorithm. 

2 Adaptive  algorithm  description 

2.1  Control  architecture 

The  model  analysis  indicates  that  the  combustion  process 
can  be  modeled  by  a feedback  interconnection  of  a second 
order  lightly  damped  system  representing  an  acoustic  mode, 


and  a large  delay  followed  by  a saturated  nonlinear  function 
representing  the  heat  release  process  [15]  [13].  If  the  phase 
shift  due  to  the  heat  release  delay  provides  a positive  feed- 
back, then  the  damping  of  the  system  will  be  reduced  and 
possibly  achieve  a negative  value.  In  the  latter  case  the  sys- 
tem dynamics  can  settle  on  a limit  cycle.  One  can  easily 
show  that  phase-shifting  control  is  an  effective  way  of  reduc- 
ing the  magnitude  of  the  limit  cycle.  The  phase  shift  can 
be  implemented  in  many  ways,  including  introducing  a delay 
in  control,  using  phase  lead  or  lag,  using  LQG  control,  and 
self-tuned  LMS  algorithms.  For  our  purpose  phase-shifting 
using  an  observer  was  selected.  The  observer-based  phase 
shifting  control  has  several  advantages:  it  is  conceptually 
simple  and  its  digital  implementation  requires  less  CPU  time 
and  offers  better  phase  resolution  that  the  delay  algorithm. 
However,  the  most  important  feature  of  the  observer  based 
phase-shifting  controller  is  its  suitability  for  tuning.  First- 
ly, the  controller  phase-shift  is  directly  available  for  tuning. 
Secondly,  an  estimate  of  the  magnitude  of  the  pressure  oscil- 
lations is  readily  available  from  the  observer  states.  Hence, 
an  extremum-seeking  algorithm  to  tune  the  controller  phase 
to  reduce  pressure  magnitude  can  be  implemented. 

There  was  a technical  challenge  in  finding  an  adaptive 
observer  that  works  for  pressure  oscillations  varying  over 
large  frequency  range  and  reliable  pressure  magnitude  de- 
tection mechanism.  This  challenge  has  been  overcome  with 
a frequency  tracking  observer  algorithm  based  on  Extended 
Kalman  Filter  described  in  [12]  that  allows  for  fast  and  reli- 
able estimation  of  the  in-phase  and  quadrature  components 
of  the  bulk  pressure  mode  over  large  range  of  frequencies. 
An  extremum-seeking  algorithm  has  been  selected  for  tun- 
ing of  the  control  phase.  One  of  the  reasons  for  choosing 
extremum-seeking  for  tuning  the  phase  is  that  it  is  possible 
to  guarantee  convergence  of  the  control  parameters  to  the 
optimal  value  and  stability  of  the  overall  system  [9,  10). 

Figure  1 shows  the  structure  of  the  control  algorithm.  The 
control  gain  was  fixed  and  only  the  control  phase  has  been 
updated.  Two  distinct  extremum-seeking  algorithms  were 
used.  The  first  one  described  in  more  details  in  this  pa- 
per relied  on  estimation  of  derivative  of  the  pressure  magni- 
tude with  respect  to  control  phase  by  introducing  a smal- 
1 sinusoidal  variation  in  the  control  phase  and  measuring 
the  response  of  the  pressure  magnitude  at  the  correspond- 
ing frequency.  The  mean  control  phase  was  incremented  by 
amount  proportional  to  the  estimated  derivative.  This  al- 
gorithm will  be  referred  to  as  a classical  extremum-seeking 
algorithm.  The  other  algorithm,  described  in  [18],  called  the 
triangular  search  algorithm,  used  three  past  sampled  average 
magnitude  values  to  determine  the  new  control  phase.  The 
idea  behind  this  algorithm  was  to  first  find  an  interval  that 
contains  the  extremum  of  a function  and  reduce  the  size  of 
this  interval  at  an  exponential  rate. 

2.2  Frequency  Tracking  Observer  Design 

Observer  based  phase-shifting  control  algorithm  that  can 
work  during  engine  acceleration  and  deceleration  transients 
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Figure  1:  Adaptive  algorithm  for  control  of  harmonic  oscil- 
lations 


requires  a frequency  tracking  algorithm.  A generalized  ver- 
sion of  a frequency  tracking  algorithm  based  on  an  Extended 
Kalman  Filter  studied  in  [12]  has  been  implemented.  Below 
we  describe  a digital  version  of  the  algorithm. 

Consider  a problem  of  estimation  of  in-phase  and  quadra- 
ture components  of  a sinusoidal  signal  y(t)  = M sin(27r  ft) 
with  unknown  frequency  /,  sampled  at  sample  time  Ts  from 
a noisy  measurement.  Th  idea  used  in  [12]  was  to  represent 
the  measurements  y(k ) = M sm(2nfkTg)  as  the  output  of  a 
nonlinear  discrete-time  system  with  the  unknown  constant 
frequency  being  one  of  the  states. 


xc(k  + 1) 

= cos(v)(h))xc(k)  — s\n(w(k))xs(k) 

(1) 

xa(k  + 1) 

= sm(w(k))xc(k)  + cos(w(k))xs(k) 

(2) 

w(k  H~  1) 

= ze(&). 

(3) 

corresponding  output  equation  is 

y(k)  = xs(k). 

(4) 

Given  measurements  y(k)  = M sm(27r  fkT8),  under  certain 
conditions  one  can  estimate  the  unknown  frequency  / (from 
the  formula  / = and  the  in-phase  and  quadrature 

components  of  the  signal  xs(k ),  xc(k)  using  an  observer. 

We  constructed  an  estimator  for  the  nonlinear  model  (1)- 
(3)  using  an  Extended  Kalman  Filter,  i.e.,  a Kalman  Filter 
designed  for  the  linearization  of  the  model  about  the  current 
estimate  (see  [1]).  In  general,  an  Extended  Kalman  Filter 
algorithm  estimates  state  x(k)  of  the  nonlinear  discrete-time 
system 

(5) 


from  the  measurement 

y(k)  = h(x{k))  + v2(k).  (6) 

ui(A;),  v2(k)  are  assumed  to  be  Gaussian,  independent,  white 
noise  signals,  with  zero  means,  and  covariance  matrices  Q, 
R.  respectively. 

The  Extended  Kalman  Filter  (to  be  sometimes  abbreviat- 
ed to  EKF)  uses  the  measurement  y(k)  to  update  the  esti- 
mate of  the  state.  The  filter  state  is  denoted  x(k)  and  the 
state  estimate  error  covariance  matrix  is  denoted  Ps(k).  The 
filter  states  are  initialized  at  zero  except  for  the  frequency 
estimate. 

Extended  Kalman  Filter  update  equations  are: 

given:  previous  filter  state  x(k  — 1),  state  estimate  er- 

ror covariance  matrix  P8(k  — 1),  measured  signal  y{k). 

1.  K(k)  = P8{k  - l)H{k  - 1 Y * (H(k  - 1 )Ps(k  - 
l)H(h  — iy  + #)-1  calculate  update  gain 

2.  yd(k)  = y(k)  — h(x(k  — 1))  filter  output  prediction 
error 

3.  £m(A:)  = x(k  — 1)  -F  K{k)y<i(k)  state  estimate  using 
prediction  error 

4.  H{k)  = , at  x = x(k  — 1)  linearization  of  out- 

put function 

5.  F(k)  = §£,,  at  x — xm(k)  linearization  of  vector 
field  / 

6.  Pm(k)  = (/  — K(k)H(k))P8{k  — 1)  state  estimate 
error  covariance  update 

7.  x(k)  ~ f{xm(h))  state  prediction 

8.  Ps(k)  = F(k)Pm(k)F(ky  H-  Q state  estimate  error 
covariance  prediction 

We  assume  R = 1 and  Q = diag{qu . . . ,<?n).  The 
choice  of  weights  </,  is  the  main  tuning  mechanism  for  the 
EKF.  While  they  were  introduced  in  the  formulation  of  the 
filtering  problem  was  as  the  variances  to  the  state  noise 
terms  they  can  be  thought  of  as  tuning  knobs  for  the  filter. 

A modified  and  extended  version  of  the  algorithm  from 
[12]  has  been  implemented  in  a digital  control  system  on  a D- 
SP  board  to  provide  a frequency  tracking  and  phase- shifting 
mechanism  for  control  of  combustion  instability.  Several  ex- 
tensions of  model  (1-2)  were  implemented.  The  most  com- 
plex one  has  8 states  and  involves  modeling  of  first  two  har- 
monics of  the  bulk  pressure  mode,  its  unknown  frequency, 
DC  (low  frequency)  component,  and  longitudinal  acoustic 
(750Hz)  mode  component. 

Figures  2 and  3 illustrate  performance  of  the  frequency 
tracking  Extended  Kalman  Filter.  An  8-second  long  data 
record  has  been  constructed  off-line  by  concatenating  four  2- 


x(k  -h  1)  = /(x(fc))  + ui(fc) 
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second  long  time  traces  of  a pressure  signal  obtained  during 
control  experiment  on  United  Technologies  Research  Center 
test  rig.  The  four  data  points  correspond  to  four  different 
control  phases.  Figure  2 shows  that  the  frequency  estimate 
from  the  EKF  algorithm  agrees  well  the  one  obtained  with 
a short-time  PSD  analysis.  From  Figure  3 we  see  that  the 
EKF  output  traces  the  pressure  signal  very  well. 


Figure  2:  Frequency  estimate  from  EKF  and 

short-time  PSD 
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Figure  3:  Pressure  signal  and  its  reconstruction 
with  EKF 

Figures  4 to  9 show  frequency  tracking  with  Extended 
Kalman  Filter  when  the  combustor  in  the  sector  rig  was 
warming  up.  Figures  4 to  6 show  the  on-line  frequency  esti- 
mate from  EKF  and  short-time  PSD.  The  frequency  of  the 
bulk  acoustic  mode  increases  because  the  speed  of  sound  in- 
creases with  temperature.  Note  that  the  frequency  tracking 
algorithm  was  tracking  the  maximum  frequency  of  the  bulk 
pressure  mode.  The  frequency  estimate  from  EKF  algorith- 
m was  smoother  than  the  one  shown  in  the  previous  section 
because  the  estimate  from  EKF  was  smoothen  with  an  addi- 
tional low-pass  filter.  Figures  7 to  9 show  pressure  signal  and 
on-line  estimates  of  the  in-phase  component  of  the  pressure 
signal  reconstructed  with  EKF,  From  comparison  of  FFTs 
one  can  see  the  band-pass  filtering  properties  of  the  EKF. 
The  in-phase  and  quadrature  components  of  the  bulk  pres- 
sure mode  have  spectral  content  around  20QHz  and  thus  the 
phase-shifted  control  signal  would  not  be  contaminated  with 
other  frequencies. 


Figures  4 to  9 illustrate  how  fast  the  pressure  oscillations 
frequency  and  magnitude  can  increase  when  the  combustor 
operating  conditions  change.  Figure  4 alone  shows  increase 
of  frequency  of  pressure  oscillation  by  30  Hz  in  9 seconds. 
Overall,  in  less  than  2 minutes  the  frequency  of  oscillations 
changed  from  about  150Hz  to  around  230Hz  and  its  mag- 
nitude increased  by  a factor  of  3.  Changes  in  the  frequen- 
cy and  magnitude  of  acoustic  oscillations  during  an  engine 
acceleration  or  deceleration  could  be  even  faster  and  more 
dramatic.  It  is  very  unlikely  that  a robust  controller  could 
cope  with  such  conditions.  The  need  for  controller  phase  and 
gain  scheduling  with  operating  conditions  is  clear.  Since  for 
the  lack  of  models  the  model-based  scheduling  of  controller 
parameters  is  impossible,  scheduling  of  controller  parameters 
using  an  adaptive  algorithm  remains  the  only  option. 

2.3  Controller  phase  tuning  using  classical 

extremum-seeking  algorithm 

The  classic  algorithm  for  finding  control  phase  9 such  that  the 
pressure  magnitude  achieves  a local  minimum  at  9*  relies  on 
introducing  a sinusoidal  variation  into  the  control  parameter 
6.  Namely,  at  time  t the  control  parameter  is  given  by 

0(t)  = 6m(t)  + Q\  sin (ojt).  (7) 

9m(t)  is  the  mean  value  of  the  control  parameter  that  is  ad- 
justed using  an  update  equation  that  is  described  later.  9\  is 
the  control  parameter  variation  magnitude,  usually  chosen  to 
be  small.  If  the  frequency  u is  small  enough,  the  magnitude 
£o(£)  approximately  follows  the  value  of  corresponding  to 
the  magnitude  of  the  limit  cycle  corresponding  to  the  control 
phase,  so  that  the  measured  magnitude  is 

A(t)  » g(Om(t)  + 6\  sin (ut))  + v(t),  (8) 

where  i/(t)  is  a random  component  of  the  pressure  magni- 
tude that  models  disturbances  driving  the  combustion  pro- 
cess, and  g(.)  is  the  function  representing  the  dependence  of 
the  magnitude  of  the  limit  cycle  corresponding  on  the  control 
phase.  Assuming  that  9\  is  small,  one  can  approximate  .4(t) 
as 

A(t)  « #.(«))  + ^ (*m(m  sin M)  + u(t).  (9) 

An  observer  has  been  used  to  extract  the  components  of  the 
magnitude  estimate  of  A(t)  at  the  frequency  w.  The  observer 
has  been  constructed  assuming  that  A(t)  is  composed  of  a 
constant  component  and  a harmonic  component  at  the  fre- 
quency a;,  i.e.,  the  magnitude  estimate  A(t)  is  an  output  of 
a linear  system  with  three  states,  one  modeling  the  constant 
component,  and  two  modeling  in-phase  and  quadrature  com- 
ponent of  a harmonic  oscillator  with  the  frequency  to.  The 
observer  equations  were 


Xi(t) 

= -cjX'z(t)  + h(A(t)  - A0(t)) 

(10) 

x2(t) 

= uxi(t)  + h(A{t)  - Aa(t)) 

(U) 

&3  (t) 

= l3(A(t)  - A„(t)) 

(12) 

A0(t) 

= x2{t)  + x3{t). 

(13) 
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with  the  gains  h,  k , h chosen  so  that  the  observer  error 
dynamics  is  stable.  (In  fact,  a discrete-time  version  of  this 
observer  was  implemented).  Now,  assuming  that  the  observ- 
er gains  are  properly  chosen,  one  can  assume  that  after  some 
transient  the  observer  states  are  approximately 


Xl(t)  s 

S (^{6,n(t))ex  eosM) 

(14) 

x2(t)  s 

a sinM) 

(15) 

x3(t)  s 

s 9(0m(t))- 

(16) 

Now,  an  estimate  of  gf(#m(£))  can  be  obtained  by  the  fol- 
lowing demodulation  of  observer  states: 

z(t)  ar2(t)  sin(wf)  + aq  (*)  cos(wt)  « ^-(0m(*))0i.  (17) 

The  control  parameter  equation  will  take  form 

4 n(t)  ~ ~az(t),  (18) 

The  parameters  to  choose  are:  the  control  parameter  vari- 
ation frequency  w,  the  observer  gains  1%,  h , h ? and  the 
update  gain  — <7. 

2.4  Controller  phase  tuning  using  triangu- 
lar search  algorithm 

The  triangular  search  algorithm  is  a self-exciting  scheme.  No 
external  probing  is  required  to  find  the  direction  of  search. 
The  purpose  of  the  contraction  coefficient  is  to  gradually  de- 
crease the  level  of  excitation  when  the  control  phase  9 is  close 
to  the  optimizing  point.  In  a noise  free  and  time  invariant 
environment,  this  is  to  guarantee  convergence  with  no  resid- 
ual error  such  as  the  case  of  a region  of  attraction.  However, 
when  noise  is  present  and  the  plant  is  time  varying,  the  ex- 
citation level  must  be  kept  at  a sufficiently  large  level  for  all 
time  in  order  to  extract  the  correct  information  about  the  di- 
rection along  which  to  move  the  control  input.  More  details 
on  the  algorithm  are  presented  in  [18]. 

3 Experimental  Results 

In  this  section  we  present  results  of  experiments  in  United 
Technologies  Research  Center  conducted  on  4 MW  Single 
Nozzle  Rig  in  August  1998  using  full-scale  engine  fuel  nozzle 
at  realistic  operating  conditions.  About  10%  of  the  net  fuel 
was  modulated  for  control  purposes  using  linear  proportional 
valve.  (For  more  details  on  the  UTRC  experiments  see  [6] 
and  [8].  ) 

Performance  specifications  for  the  adaptive  algorithm  have 
been  defined  for  algorithm  initialization  transients  and  en- 
gine acceleration  transients.  When  initialized  with  a phase 
corresponding  to  amplification  of  oscillations,  the  algorithms 
should  quickly  produce  and  maintain  phases  corresponding 
to  suppression  of  the  oscillations.  During  the  engine  accel- 
eration transients  the  algorithms  should  be  able  to  suppress 
oscillations  relative  to  uncontrolled  levels. 


To  test  the  transient  performance  of  the  adaptive  algorith- 
m initialization  transients  have  been  introduced.  The  initial 
control  phase  varied  significantly  from  the  optimal  one.  The 
algorithms  behaved  very  well  at  higher  power  condition  (s- 
iriall  pressure  oscillations  and  medium  level  of  broad-band 
disturbance)  and  reasonably  well  at  lower  power  condition- 
s (large  pressure  oscillations  and  large  level  of  broad-band 
disturbance) . Once  reaching  a neighborhood  of  the  optimal 
value,  the  control  phase  usually  stayed  in  a reasonably  small 
neighborhood  of  that  value,  rarely  produced  control  phases 
corresponding  to  level  higher  than  uncontrolled  levels,  and 
always  provided  better  average  pressure  oscillations  levels 
than  uncontrolled  levels. 

The  triangular  search  algorithm  tracked  the  minimum  val- 
ue well  in  a transient  from  higher  power  to  lower  power  con- 
dition. This  was  a slow  transient  but  the  frequency  changed 
by  about  40Hz  and  the  pressure  magnitude  and  noise  levels 
were  changing  dramatically.  Classical  algorithm  with  a fixed 
gain  did  not  work  during  the  transients,  as  the  controller  gain 
that  was  optimal  at  the  higher  power  level  was  destabilizing 
the  system  at  the  lower  power  level.  Gain  adaptation  would 
be  needed  for  the  classical  algorithm  to  work  at  both  pow- 
er levels.  The  frequency  and  magnitude  tracking  mechanism 
worked  very  well  during  the  transients. 

The  dependence  of  the  mean  pressure  magnitude  and  fre- 
quency of  the  corresponding  mode  on  the  control  phase  has 
been  determined  experimentally  at  two  power  conditions,  so 
that  the  optimal  control  phase  was  known  a priori.  This 
information  allowed  to  check  performance  of  the  extremum- 
seeking algorithms.  Initialization  transients  have  been  in- 
troduced, where  the  initial  control  phase  varied  significantly 
from  the  optimal  one.  Figures  10  to  31  show  typical  time 
histories  of  control  phase  and  pressure  magnitude  (instan- 
taneous and  mean)  estimate  at  two  power  conditions.  The 
open-loop  and  optimal  levels  of  the  pressure  as  well  as  the 
optimal  control  phase  are  indicated  by  horizontal  lines  in  the 
corresponding  plots,  oscillations 

The  experimental  results  are  presented  in  Appendix. 

For  more  details  on  experiments  we  refer  to  [5].  In  general 
both  phase  tuning  algorithms  behaved  very  well  at  higher 
power  condition  (medium  noise  and  small  pressure  oscilla- 
tios)  and  reasonably  well  at  lower  power  conditions  (large 
noise  and  pressure  oscillations).  Once  reaching  a neighbor- 
hood of  the  optimal  value,  the  control  phase  usually  stayed 
in  a reasonably  small  neighborhood  of  that  value,  rarely  pro- 
duced control  phases  corresponding  to  level  higher  than  un- 
controlled levels,  and  always  provided  better  average  pressure 
oscillations  levels  than  uncontrolled  levels. 

It  has  been  concluded  that  the  major  factor  affecting  the 
performance  of  the  extremum-seeking  schemes  is  the  "noise” 
present  in  the  pressure  magnitude.  More  precisely,  the  main 
trouble  with  implementation  of  extremum-seeking  algorithm 
is  that  the  pressure  oscillations  do  not  have  a consistent  in- 
stantaneous magnitude.  In  fact,  the  magnitude  varies  wildly. 
However,  changing  the  control  phase  affects  the  mean  value 
of  the  pressure  and  after  some  averaging  time  it  is  possible 
to  determine  the  mean  pressure  magnitude  as  a function  of 
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control  phase  with  one  minimum  and  one  maximum  (modu- 
lo 360  degrees).  The  fluctuating  component  of  the  pressure 
magnitude  for  a fixed  phase  is  what  we  call  noise.  The  noise 
can  be  attributed  to  an  efect  of  turbulent  flow  in  the  com- 
bustor on  the  acoustic  mode  either  directly,  or  via  chemical 
reaction. 

Since  magnitude  of  the  noise  is  a crucial  factor  that  deter- 
mines performance  of  the  timing  mechanism,  it  is  necessary 
to  include  a representative  noise  model  in  the  engine  dynam- 
ics model,  if  the  latter  is  used  to  study  performance  of  the 
tuning  algorithms. 

3.1  Simulated  engine  acceleration  transient 

A nominal  gas  turbine  engine  transient  from  50%  to  100% 
power  conditions  lasts  a fraction  of  a minute.  The  changes 
in  operating  conditions  appearing  during  engine  acceleration 
and  deceleration  are  likely  to  resemble  the  transients  between 
different  power  levels  on  the  single  nozzle  rig.  Engine  tran- 
sient lasting  fraction  of  a minute  could  not  be  simulated 
in  experiments  on  the  single  nozzle  rig.  Instead,  we  test- 
ed the  extremum-seeking  algorithms  in  a simulated  engine 
30-second  transient,  where  the  measured  magnitude  versus 
control  phase  functions  and  lower  and  higher  power  were  in- 
terpolated. 

To  simulate  the  transients  the  form  of  the  averaged  pres- 
sure magnitude  model  was  assumed  in  accordance  with  the 
predictions  of  analysis  of  the  phvsics-based  reduced  order 
model  and  the  unknown  parameters  of  the  averaged  mod- 
el were  determined  from  experiments  (experiments  are  de- 
scribed in  [5]).  The  identified  model  was 

x0(t)  = -a(0(t))  {xo(t)  - g(0(t)))  (19) 

A{t)  = xo  (t)  + i/(t),  (20) 

where  0(t)  is  the  controller  phase,  xo  (t)  is  the  mean  pressure 
magnitude,  A(t)  is  the  instantaneous  pressure  magnitude, 
and  u(t)  is  a random  noise.  The  functions  <?(.),  a(.),  and  the 
statistics  of  (colored)  noise  v(t)  were  identified  from  exper- 
iment. To  simplify  analysis  we  do  not  include  the  frequen- 
cy tracking  algorithm  in  the  model.  This  is  justified  under 
the  assumption  that  the  frequency  changes  faster  than  the 
magnitude  and  that  the  frequency  tracking  algorithm  is  fast 
and  stable.  In  other  words,  we  assume  that  the  frequency 
dynamics  and  the  tracking  algorithm  contribute  to  the  fast 
stable  dynamics  whose  time  scale  is  separated  from  the  slow 
magnitude  dynamics. 

The  triangular  search  algorithm  performed  well  in  the  sim- 
ulated transient  (see  Figure  11).  It  was  determined  that  in 
order  to  work  in  a simulated  transient  from  lower  to  higher 
power  conditions,  the  classical  algorithm  would  have  to  be 
modified  to  allow  for  adaptive  gain  change  (by  a factor  of 
five).  One  fixed  gain  would  not  work  at  both  power  condi- 
tions. 

Ability  of  an  extremum-seeking  algorithm  to  track  the  fast 
changing  optimal  control  phase  shift  during  the  engine  accel- 
eration and  deceleration  transient  conditions  in  the  present 
of  disturbance  driving  the  system  has  been  demonstrated  in 


a simulation.  However,  there  is  need  to  study  stability,  ro- 
bustness, and  performance  of  the  algorithms  using  analytical 
tools.  Traditional  methods  based  on  time-scale  separation 
are  not  sufficient  as  the  time  scale  of  of  change  of  operating 
conditions  is  not  well  separated  from  the  time  scale  of  the 
transients  in  the  dynamics.  Preliminary  analysis  of  perfor- 
mance limiations  of  extremum-seeking  algorithms  that  does 
not  expoits  time-scale  separation  has  been  presented  in  [10]. 

Conclusions 

We  presented  results  of  experiments  with  two  distinct 
extremum-seeking  adaptive  algorithms  for  control  of  com- 
bustion instability  suitable  for  reduction  of  acoustic  pressure 
oscillations  in  gas  turbine  over  large  range  of  operating  con- 
ditions. The  algorithms  consists  of  a frequency  tracking  Ex- 
tended Kalman  Filter  to  determine  the  in-phase  component, 
the  quadrature  component,  and  the  magnitude  of  the  acous- 
tic mode  of  interest,  and  a phase  shifting  controller  with 
the  controller  phase  tuned  using  an  extremum-seeking  al- 
gorithms. The  algorithms  have  been  tested  in  combustion 
experiments  involving  full-scale  engine  hardware  and  during 
simulated  fast  engine  transients. 
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Figure  4:  Frequency  estimate 

from  EKF  and  short-time  PSD 
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Figure  7;  Pressure  signal  and  its 
in-phase  component  reeoritr  acted 
with  EKF:  time  traces  and  FFT 
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Figure  5:  Frequency  estimate 

from  EKF  and  short-time  PSD 
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Figure  8:  Pressure  signal  and  its 
in-phase  component  recontracted 
with  EKF:  time  traces  and  FFT 
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Figure  6:  Frequency  estimate 

from  EKF  and  short-time  PSD 


Figure  9:  Pressure  signal  and  its 
In-phase  component  recontracted 
with  EKF:  time  traces  and  FFT 
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Small  periodic  variation  is  ksirodueed  in  control  phase  to  estimate 
derivative  of  pressure  magnitude  with  respect  to  control  phase 
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Figure  10:  Classical  extremum-seeking  algorithm.  Typical 
initialization  transient  for  higher  power.  Control  phase  and 
pressure  magnitude  as  functions  of  time.  Horizontal  lines 
show  optimal  and  close  loop  levels  of  mean  pressure  magni- 
tude. 
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Figure  12:  Controller  phase  tun- 
ing using  classical  algorithm.  Ini- 
tialization transient  for  higher 
power.  / = 1 Hz,  9\  = 15, 
k = 1000,  M0)  = 170.  Control 
phase  and  pressure  magnitude  as 
functions  of  time. 
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Figure  1 1 : Perfomance  of  triangular  extremum-seeking  algo- 
rithm during  a simulated  engine  transients 


Figure  13:  Controller  phase  tun- 
ing using  classical  algorithm.  Ini- 
tialization transient  for  higher 
power.  / = 1 Hz,  9\  = 15, 
k = 1000,  M0)  = -10.  Control 
phase  and  pressure  magnitude  as 
functions  of  time. 
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Figure  14:  Controller  phase  tun- 
ing using  classical  algorithm.  Ini- 
tialization transient  for  higher 
power.  / = 1 Hz,  9i  — 15,  k ~ 
1000,  0TO(O)  = 60.  Control  phase 
and  pressure  magnitude  as  func- 
tions of  time. 
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Figure  15:  Controller  phase  tun- 
ing using  classical  algorithm-  Ini- 
tialization transient  for  higher 
power.  / = 1 Hz,  9\  = 15, 
k = 1000,  #m(0)  = 170  . Pressure 
magnitude  as  function  of  control 
phase. 
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Figure  18:  Controller  phase  tun- 
ing using  classical  algorithm.  Ini- 
tialization transient  for  lower 
power.  / = 1 Hz,  9\  = 15, 
k = 150,  0m(O)  = —135.  Control 
phase  and  pressure  magnitude  as 
functions  of  time. 
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Figure  16:  Controller  phase  tun- 
ing using  classical  algorithm.  Ini- 
tialization transient  for  higher 
power.  / = 1 Hz,  9\  ~ 15, 
k = 1000,  0m{ 0)  = -10.  Pressure 
magnitude  as  function  of  control 
phase. 


Figure  19:  Controller  phase  tun- 
ing using  classical  algorithm.  Ini- 
tialization transient  for  lower 
power.  / = 1 Hz,  9 ] = 15,  k = 
150,  0m(O)  = 45.  Control  phase 
and  pressure  magnitude  as  func- 
tions of  time. 


Figure  17:  Controller  phase  tun- 
ing using  classical  algorithm.  Ini- 
tialization transient  for  higher 
power.  / = 1 Hz,  9 \ = 15, 
k = 1000,  9m( 0)  = 60.  Pressure 
magnitude  as  function  of  control 
phase. 


Figure  20:  Controller  phase  tun- 
ing using  classical  algorithm.  Ini- 
tialization transient  for  lower 
power.  / = lHz:  9i  = 15,  k — 
150,  0m(O)  = 115.  Control  phase 
and  pressure  magnitude  as  func- 
tions of  time. 
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Figure  21:  Controller  phase  tun- 
ing using  classical  algorithm.  Ini- 
tialization transient  for  lower 
power.  / = 1 Hz,  9\  = 15,  k = 
150,  MO)  = -135  . Pressure 
magnitude  as  function  of  control 
phase. 


Figure  22:  Controller  phase  tun- 
ing using  classical  algorithm.  Ini- 
tialization transient  for  lower 
power.  / = 1Hz,  &i  = 15, 
k = 150,  0m(O)  = 45.  Pressure 
magnitude  as  function  of  control 
phase. 
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Figure  23:  Controller  phase  tun- 
ing using  classical  algorithm.  Ini- 
tialization transient  for  lower 
power.  / = IFjTz,  6i  = 15, 
k — 150,  6m{ 0)  = 115.  Pressure 
magnitude  as  function  of  control 
phase. 


Figure  24:  Controller  phase  tun- 
ing using  triangular  search  algo- 
rithm. Initialization  transient  for 
higher  power.  0m(O)  = 170.  Con- 
trol phase  and  pressure  magni- 
tude as  functions  of  time. 
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Figure  25:  Controller  phase  tun- 
ing using  triangular  search  algo- 
rithm. Initialization  transient  for 
higher  power.  0m(0)  = -10. 
Control  phase  and  pressure  mag- 
nitude as  functions  of  time. 
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Figure  26:  Controller  phase  tun- 
ing using  triangular  search  algo- 
rithm. Initialization  transient  for 
higher  powder.  0m(O)  = 60.  Con- 
trol phase  and  pressure  magni- 
tude as  functions  of  time. 
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Figure  27:  Controller  phase  tun- 
ing using  triangular  search  algo- 
rithm. Initialization  transient  for 
higher  power.  i9m(0)  = 170  . 
Pressure  magnitude  as  function 
of  control  phase. 
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Figure  28:  Controller  phase  tun- 
ing using  triangular  search  algo- 
rithm. Initialization  transient  for 
higher  power,  0m  (0)  = -10. 
Pressure  magnitude  as  function 
of  control  phase. 
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Figure  30:  Controller  phase  tun- 
ing using  triangular  search  algo- 
rithm. Initialization  transient  for 
lower  power.  0m(O)  = -135. 
Control  phase  and  pressure  mag- 
nitude as  functions  of  time. 
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Figure  31:  Controller  phase  tun- 
ing using  triangular  search  algo- 
rithm. Initialization  transient  for 
lower  power.  0m(0)  = 45.  Con- 
trol phase  and  pressure  magni- 
tude as  functions  of  time. 


Figure  29:  Controller  phase  tun- 
ing using  triangular  search  al- 
gorithm. Initialization  transient 
for  higher  power.  0m(O)  = 60. 
Pressure  magnitude  as  function 
of  control  phase. 
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Figure  32;  Controller  phase  tun- 
ing using  triangular  search  algo- 
rithm. Initialization  transient  for 
lower  power.  . Pressure  magni- 
tude as  function  of  control  phase. 


Figure  33:  Controller  phase  tun- 
ing using  triangular  search  al- 
gorithm. Initialization  transien- 
t for  lower  power.  0m(0)  = 45. 
Pressure  magnitude  as  function 
of  control  phase. 


PAPER  -31,  A.  Banaszuk 

Question  (F.  E.  C.  Culick,  USA) 

Is  there  (i.e.,  can  you  give)  a simple  reason  why  your  adaptive  control  algorithm! s) 
worked  better  for  high  power  than  for  lower  power  conditions? 

Reply 

It  appears  that  the  high  power  condition  was  characterized  by  a linear  model  and  the  flat 
bottomed  behavior  of  the  phase-magnitude  map  was  caused  by  peak-splitting  (inability  to 
improve  performance  beyond  a certain  point).  We  speculate  that  the  “notchy”  shape  of 
the  phase-magnitude  function  is  caused  by  the  system  going  in  and  out  of  a limit  cycle. 
The  limit  cycling  behavior  was  caused  by  a subcritical  bifurcation.  This  increases  sensi- 
tivity of  the  magnitude  to  the  phase  change  around  the  minimum. 
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